Data-driven identification of predictive risk biomarkers for subgroups of osteoarthritis using interpretable machine learning.

Osteoarthritis (OA) is increasing in prevalence and has a severe impact on patients' lives. However, our understanding of biomarkers driving OA risk remains limited. We developed a model predicting the five-year risk of OA diagnosis, integrating retrospective clinical, lifestyle and biomarker data from the UK Biobank (19,120 patients with OA, ROC-AUC: 0.72, 95%CI (0.71-0.73)). Higher age, BMI and prescription of non-steroidal anti-inflammatory drugs contributed most to increased OA risk prediction ahead of diagnosis. We identified 14 subgroups of OA risk profiles. These subgroups were validated in an independent set of patients evaluating the 11-year OA risk, with 88% of patients being uniquely assigned to one of the 14 subgroups. Individual OA risk profiles were characterised by personalised biomarkers. Omics integration demonstrated the predictive importance of key OA genes and pathways (e.g., GDF5 and TGF-ß signalling) and OA-specific biomarkers (e.g., CRTAC1 and COL9A1). In summary, this work identifies opportunities for personalised OA prevention and insights into its underlying pathogenesis.

Longitudinal analysis of serum-derived extracellular vesicle RNA to monitor dacomitinib treatment response in EGFR-amplified recurrent glioblastoma patients.

Background: Glioblastoma (GBM) is a highly aggressive and invasive brain tumor associated with high patient mortality. A large fraction of GBM tumors have been identified as epidermal growth factor receptor (EGFR) amplified and ~50% also are EGFRvIII mutant positive. In a previously reported multicenter phase II study, we have described the response of recurrent GBM (rGBM) patients to dacomitinib, an EGFR tyrosine kinase inhibitor (TKI). As a continuation of that report, we leverage the tumor cargo-encapsulating extracellular vesicles (EVs) and explore their genetic composition as carriers of tumor biomarker. Methods: Serum samples were longitudinally collected from EGFR-amplified rGBM patients who clinically benefitted from dacomitinib therapy (responders) and those who did not (nonresponders), as well as from a healthy cohort of individuals. The serum EV transcriptome was evaluated to map the RNA biotype distribution and distinguish GBM disease. Results: Using long RNA sequencing, we show enriched detection of over 10 000 coding RNAs from serum EVs. The EV transcriptome yielded a unique signature that facilitates differentiation of GBM patients from healthy donors. Further analysis revealed genetic enrichment that enables stratification of responders from nonresponders prior to dacomitinib treatment as well as following administration. Conclusion: This study demonstrates that genetic composition analysis of serum EVs may aid in therapeutic stratification to identify patients with dacomitinib-responsive GBM.

Hepatocyte mARC1 promotes fatty liver disease.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) has a prevalence of ∼25% worldwide, with significant public health consequences yet few effective treatments. Human genetics can help elucidate novel biology and identify targets for new therapeutics. Genetic variants in mitochondrial amidoxime-reducing component 1 (MTARC1) have been associated with NAFLD and liver-related mortality; however, its pathophysiological role and the cell type(s) mediating these effects remain unclear. We aimed to investigate how MTARC1 exerts its effects on NAFLD by integrating human genetics with in vitro and in vivo studies of mARC1 knockdown. Methods: Analyses including multi-trait colocalisation and Mendelian randomisation were used to assess the genetic associations of MTARC1. In addition, we established an in vitro long-term primary human hepatocyte model with metabolic readouts and used the Gubra Amylin NASH (GAN)-diet non-alcoholic steatohepatitis mouse model treated with hepatocyte-specific N-acetylgalactosamine (GalNAc)-siRNA to understand the in vivo impacts of MTARC1. Results: We showed that genetic variants within the MTARC1 locus are associated with liver enzymes, liver fat, plasma lipids, and body composition, and these associations are attributable to the same causal variant (p.A165T, rs2642438 G>A), suggesting a shared mechanism. We demonstrated that increased MTARC1 mRNA had an adverse effect on these traits using Mendelian randomisation, implying therapeutic inhibition of mARC1 could be beneficial. In vitro mARC1 knockdown decreased lipid accumulation and increased triglyceride secretion, and in vivo GalNAc-siRNA-mediated knockdown of mARC1 lowered hepatic but increased plasma triglycerides. We found alterations in pathways regulating lipid metabolism and decreased secretion of 3-hydroxybutyrate upon mARC1 knockdown in vitro and in vivo. Conclusions: Collectively, our findings from human genetics, and in vitro and in vivo hepatocyte-specific mARC1 knockdown support the potential efficacy of hepatocyte-specific targeting of mARC1 for treatment of NAFLD. Impact and implications: We report that genetically predicted increases in MTARC1 mRNA associate with poor liver health. Furthermore, knockdown of mARC1 reduces hepatic steatosis in primary human hepatocytes and a murine NASH model. Together, these findings further underscore the therapeutic potential of targeting hepatocyte MTARC1 for NAFLD.

Exosome based analysis for Space Associated Neuro-Ocular Syndrome and health risks in space exploration.

Molecular profiling to characterize the effects of environmental exposures is important from the human health and performance as well as the occupational medicine perspective in space exploration. We have developed a novel exosome-based platform that allows profiling of biological processes in the body from a variety of body fluids. The technology is suitable for diagnostic applications as well as studying the pathophysiology of the Space Associated Neuro-Ocular Syndrome in astronauts and monitoring patients with chronically impaired cerebrospinal fluid drainage or elevated intracranial pressure. In this proof-of-concept, we demonstrate that: (a) exosomes from different biofluids contain a specific population of RNA transcripts; (b) urine collection hardware aboard the ISS is compatible with exosome gene expression technology; (c) cDNA libraries from exosomal RNA can be stored in dry form and at room temperature, representing an interesting option for the creation of longitudinal molecular catalogs that can be stored as a repository for retrospective analysis.

Phase 2 of extracellular RNA communication consortium charts next-generation approaches for extracellular RNA research.

The extracellular RNA communication consortium (ERCC) is an NIH-funded program aiming to promote the development of new technologies, resources, and knowledge about exRNAs and their carriers. After Phase 1 (2013-2018), Phase 2 of the program (ERCC2, 2019-2023) aims to fill critical gaps in knowledge and technology to enable rigorous and reproducible methods for separation and characterization of both bulk populations of exRNA carriers and single EVs. ERCC2 investigators are also developing new bioinformatic pipelines to promote data integration through the exRNA atlas database. ERCC2 has established several Working Groups (Resource Sharing, Reagent Development, Data Analysis and Coordination, Technology Development, nomenclature, and Scientific Outreach) to promote collaboration between ERCC2 members and the broader scientific community. We expect that ERCC2's current and future achievements will significantly improve our understanding of exRNA biology and the development of accurate and efficient exRNA-based diagnostic, prognostic, and theranostic biomarker assays.

Technical Performance Evaluation of Olink Proximity Extension Assay for Blood-Based Biomarker Discovery in Longitudinal Studies of Alzheimer's Disease.

The core Alzheimer's disease (AD) cerebrospinal fluid (CSF) biomarkers; amyloid-ß (Aß), total tau (t-tau), and phosphorylated tau (p-tau181), are strong indicators of the presence of AD pathology, but do not correlate well with disease progression, and can be difficult to implement in longitudinal studies where repeat biofluid sampling is required. As a result, blood-based biomarkers are increasingly being sought as alternatives. In this study, we aimed to evaluate a promising blood biomarker discovery technology, Olink Proximity Extension Assays for technical reproducibility characteristics in order to highlight the advantages and disadvantages of using this technology in biomarker discovery in AD. We evaluated the performance of five Olink Proteomic multiplex proximity extension assays (PEA) in plasma samples. Three technical control samples included on each plate allowed calculation of technical variability. Biotemporal stability was measured in three sequential annual samples from 54 individuals with and without AD. Coefficients of variation (CVs), analysis of variance (ANOVA), and variance component analyses were used to quantify technical and individual variation over time. We show that overall, Olink assays are technically robust, with the largest experimental variation stemming from biological differences between individuals for most analytes. As a powerful illustration of one of the potential pitfalls of using a multi-plexed technology for discovery, we performed power calculations using the baseline samples to demonstrate the size of study required to overcome the need for multiple test correction with this technology. We show that the power of moderate effect size proteins was strongly reduced, and as a result investigators should strongly consider pooling resources to perform larger studies using this multiplexed technique where possible.

Plasma Proteomics of COVID-19-Associated Cardiovascular Complications: Implications for Pathophysiology and Therapeutics.

To gain insights into the mechanisms driving cardiovascular complications in COVID-19, we performed a case-control plasma proteomics study in COVID-19 patients. Our results identify the senescence-associated secretory phenotype, a marker of biological aging, as the dominant process associated with disease severity and cardiac involvement. FSTL3, an indicator of senescence-promoting Activin/TGFß signaling, and ADAMTS13, the von Willebrand Factor-cleaving protease whose loss-of-function causes microvascular thrombosis, were among the proteins most strongly associated with myocardial stress and injury. Findings were validated in a larger COVID-19 patient cohort and the hamster COVID-19 model, providing new insights into the pathophysiology of COVID-19 cardiovascular complications with therapeutic implications.

Proteomic characterization of post-mortem human brain tissue following ultracentrifugation-based subcellular fractionation.

Proteomic characterization of human brain tissue is increasingly utilized to identify potential novel biomarkers and drug targets for a variety of neurological diseases. In whole-tissue studies, results may be driven by changes in the proportion of the largest and most abundant organelles or tissue cell-type composition. Spatial proteomics approaches enhance our knowledge of disease mechanisms and changing signalling pathways at the subcellular level by taking into account the importance of cellular localization, which critically influences protein function. Density gradient-based ultracentrifugation methods allow for subcellular fractionation and have been utilized in cell lines, mouse and human brain tissue to quantify thousands of proteins in specific enriched organelles such as the pre- and post-synapse. Serial ultracentrifugation methods allow for the analysis of multiple cellular organelles from the same biological sample, and to our knowledge have not been previously applied to frozen post-mortem human brain tissue. The use of frozen human tissue for tissue fractionation faces two major challenges, the post-mortem interval, during which proteins may leach from their usual location into the cytosol, and freezing, which results in membrane breakdown. Despite these challenges, in this proof-of-concept study, we show that the majority of proteins segregate reproducibly into crude density-based centrifugation fractions, that the fractions are enriched for the appropriate organellar markers and that significant differences in protein localization can be observed between tissue from individuals with Alzheimer's disease and control individuals.

Exercise hormone irisin is a critical regulator of cognitive function.

Identifying secreted mediators that drive the cognitive benefits of exercise holds great promise for the treatment of cognitive decline in ageing or Alzheimer's disease (AD). Here, we show that irisin, the cleaved and circulating form of the exercise-induced membrane protein FNDC5, is sufficient to confer the benefits of exercise on cognitive function. Genetic deletion of Fndc5/irisin (global Fndc5 knock-out (KO) mice; F5KO) impairs cognitive function in exercise, ageing and AD. Diminished pattern separation in F5KO mice can be rescued by delivering irisin directly into the dentate gyrus, suggesting that irisin is the active moiety. In F5KO mice, adult-born neurons in the dentate gyrus are morphologically, transcriptionally and functionally abnormal. Importantly, elevation of circulating irisin levels by peripheral delivery of irisin via adeno-associated viral overexpression in the liver results in enrichment of central irisin and is sufficient to improve both the cognitive deficit and neuropathology in AD mouse models. Irisin is a crucial regulator of the cognitive benefits of exercise and is a potential therapeutic agent for treating cognitive disorders including AD.

Evaluation of serological lateral flow assays for severe acute respiratory syndrome coronavirus-2.

BACKGROUND: COVID-19 has resulted in significant morbidity and mortality worldwide. Lateral flow assays can detect anti-Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) antibodies to monitor transmission. However, standardized evaluation of their accuracy and tools to aid in interpreting results are needed. METHODS: We evaluated 20 IgG and IgM assays selected from available tests in April 2020. We evaluated the assays' performance using 56 pre-pandemic negative and 56 SARS-CoV-2-positive plasma samples, collected 10-40 days after symptom onset, confirmed by a molecular test and analyzed by an ultra-sensitive immunoassay. Finally, we developed a user-friendly web app to extrapolate the positive predictive values based on their accuracy and local prevalence. RESULTS: Combined IgG + IgM sensitivities ranged from 33.9 to 94.6%, while combined specificities ranged from 92.6 to 100%. The highest sensitivities were detected in Lumiquick for IgG (98.2%), BioHit for both IgM (96.4%), and combined IgG + IgM sensitivity (94.6%). Furthermore, 11 LFAs and 8 LFAs showed perfect specificity for IgG and IgM, respectively, with 15 LFAs showing perfect combined IgG + IgM specificity. Lumiquick had the lowest estimated limit-of-detection (LOD) (0.1 µg/mL), followed by a similar LOD of 1.5 µg/mL for CareHealth, Cellex, KHB, and Vivachek. CONCLUSION: We provide a public resource of the accuracy of select lateral flow assays with potential for home testing. The cost-effectiveness, scalable manufacturing process, and suitability for self-testing makes LFAs an attractive option for monitoring disease prevalence and assessing vaccine responsiveness. Our web tool provides an easy-to-use interface to demonstrate the impact of prevalence and test accuracy on the positive predictive values.

Synaptic proteins associated with cognitive performance and neuropathology in older humans revealed by multiplexed fractionated proteomics.

Alzheimer's disease (AD) is defined by the presence of abundant amyloid-ß (Aß) and tau neuropathology. While this neuropathology is necessary for AD diagnosis, it is not sufficient for causing cognitive impairment. Up to one third of community dwelling older adults harbor intermediate to high levels of AD neuropathology at death yet demonstrate no significant cognitive impairment. Conversely, there are individuals who exhibit dementia with no gross explanatory neuropathology. In prior studies, synapse loss correlated with cognitive impairment. To understand how synaptic composition changes in relation to neuropathology and cognition, multiplexed liquid chromatography mass-spectrometry was used to quantify enriched synaptic proteins from the parietal association cortex of 100 subjects with contrasting levels of AD pathology and cognitive performance. 123 unique proteins were significantly associated with diagnostic category. Functional analysis showed enrichment of serotonin release and oxidative phosphorylation categories in normal (cognitively unimpaired, low neuropathology) and "resilient" (unimpaired despite AD pathology) individuals. In contrast, frail individuals, (low pathology, impaired cognition) showed a metabolic shift towards glycolysis and increased presence of proteasome subunits.

Whole-genome sequencing of phenotypically distinct inflammatory breast cancers reveals similar genomic alterations to non-inflammatory breast cancers.

BACKGROUND: Inflammatory breast cancer (IBC) has a highly invasive and metastatic phenotype. However, little is known about its genetic drivers. To address this, we report the largest cohort of whole-genome sequencing (WGS) of IBC cases. METHODS: We performed WGS of 20 IBC samples and paired normal blood DNA to identify genomic alterations. For comparison, we used 23 matched non-IBC samples from the Cancer Genome Atlas Program (TCGA). We also validated our findings using WGS data from the International Cancer Genome Consortium (ICGC) and the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We examined a wide selection of genomic features to search for differences between IBC and conventional breast cancer. These include (i) somatic and germline single-nucleotide variants (SNVs), in both coding and non-coding regions; (ii) the mutational signature and the clonal architecture derived from these SNVs; (iii) copy number and structural variants (CNVs and SVs); and (iv) non-human sequence in the tumors (i.e., exogenous sequences of bacterial origin). RESULTS: Overall, IBC has similar genomic characteristics to non-IBC, including specific alterations, overall mutational load and signature, and tumor heterogeneity. In particular, we observed similar mutation frequencies between IBC and non-IBC, for each gene and most cancer-related pathways. Moreover, we found no exogenous sequences of infectious agents specific to IBC samples. Even though we could not find any strongly statistically distinguishing genomic features between the two groups, we did find some suggestive differences in IBC: (i) The MAST2 gene was more frequently mutated (20% IBC vs. 0% non-IBC). (ii) The TGF ß pathway was more frequently disrupted by germline SNVs (50% vs. 13%). (iii) Different copy number profiles were observed in several genomic regions harboring cancer genes. (iv) Complex SVs were more frequent. (v) The clonal architecture was simpler, suggesting more homogenous tumor-evolutionary lineages. CONCLUSIONS: Whole-genome sequencing of IBC manifests a similar genomic architecture to non-IBC. We found no unique genomic alterations shared in just IBCs; however, subtle genomic differences were observed including germline alterations in TGFß pathway genes and somatic mutations in the MAST2 kinase that could represent potential therapeutic targets.

The Dynamic Platelet Transcriptome in Obesity and Weight Loss.

OBJECTIVE: Adiposity is associated with oxidative stress, inflammation, and glucose intolerance. Previous data suggest that platelet gene expression is associated with key cardiometabolic phenotypes, including body mass index but stable in healthy individuals over time. However, modulation of gene expression in platelets in response to metabolic shifts (eg, weight reduction) is unknown and may be important to defining mechanism. Approach and Results: Platelet RNA sequencing and aggregation were performed from 21 individuals with massive weight loss (>45 kg) following bariatric surgery. Based on RNA sequencing data, we measured the expression of 67 genes from isolated platelet RNA using high-throughput quantitative reverse transcription quantitative PCR in 1864 FHS (Framingham Heart Study) participants. Many transcripts not previously studied in platelets were differentially expressed with bariatric surgical weight loss, appeared specific to platelets (eg, not differentially expressed in leukocytes), and were enriched for a nonalcoholic fatty liver disease pathway. Platelet aggregation studies did not detect alteration in platelet function after significant weight loss. Linear regression models demonstrated several platelet genes modestly associated with cross-sectional cardiometabolic phenotypes, including body mass index. There were no associations between studied transcripts and incident diabetes or cardiovascular end points. CONCLUSIONS: In summary, while there is no change in platelet aggregation function after significant weight loss, the human platelet experiences a dramatic transcriptional shift that implicates pathways potentially relevant to improved cardiometabolic risk postweight loss (eg, nonalcoholic fatty liver disease). Further studies are needed to determine the mechanistic importance of these observations.

Open Problems in Extracellular RNA Data Analysis: Insights From an ERCC Online Workshop.

We now know RNA can survive the harsh environment of biofluids when encapsulated in vesicles or by associating with lipoproteins or RNA binding proteins. These extracellular RNA (exRNA) play a role in intercellular signaling, serve as biomarkers of disease, and form the basis of new strategies for disease treatment. The Extracellular RNA Communication Consortium (ERCC) hosted a two-day online workshop (April 19-20, 2021) on the unique challenges of exRNA data analysis. The goal was to foster an open dialog about best practices and discuss open problems in the field, focusing initially on small exRNA sequencing data. Video recordings of workshop presentations and discussions are available (https://exRNA.org/exRNAdata2021-videos/). There were three target audiences: experimentalists who generate exRNA sequencing data, computational and data scientists who work with those groups to analyze their data, and experimental and data scientists new to the field. Here we summarize issues explored during the workshop, including progress on an effort to develop an exRNA data analysis challenge to engage the community in solving some of these open problems.

Exploring Predictors of Response to Dacomitinib in EGFR-Amplified Recurrent Glioblastoma.

PURPOSE: Despite the high frequency of EGFR genetic alterations in glioblastoma (GBM), EGFR-targeted therapies have not had success in this disease. To improve the likelihood of efficacy, we targeted adult patients with recurrent GBM enriched for EGFR gene amplification, which occurs in approximately half of GBM, with dacomitinib, a second-generation, irreversible epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor that penetrates the blood-brain barrier, in a multicenter phase II trial. PATIENTS AND METHODS: We retrospectively explored whether previously described EGFR extracellular domain (ECD)-sensitizing mutations in the context of EGFR gene amplification could predict response to dacomitinib, and in a predefined subset of patients, we measured post-treatment intratumoral dacomitinib levels to verify tumor penetration. RESULTS: We found that dacomitinib effectively penetrates contrast-enhancing GBM tumors. Among all 56 treated patients, 8 (14.3%) had a clinical benefit as defined by a duration of treatment of at least 6 months, of whom 5 (8.9%) remained progression free for at least 1 year. Presence of EGFRvIII or EGFR ECD missense mutation was not associated with clinical benefit. We evaluated the pretreatment transcriptome in circulating extracellular vesicles (EVs) by RNA sequencing in a subset of patients and identified a signature that distinguished patients who had durable benefit versus those with rapid progression. CONCLUSION: While dacomitinib was not effective in most patients with EGFR-amplified GBM, a subset experienced a durable, clinically meaningful benefit. Moreover, EGFRvIII and EGFR ECD mutation status in archival tumors did not predict clinical benefit. RNA signatures in circulating EVs may warrant investigation as biomarkers of dacomitinib efficacy in GBM.

exceRpt: A Comprehensive Analytic Platform for Extracellular RNA Profiling.

Small RNA sequencing has been widely adopted to study the diversity of extracellular RNAs (exRNAs) in biofluids; however, the analysis of exRNA samples can be challenging: they are vulnerable to contamination and artifacts from different isolation techniques, present in lower concentrations than cellular RNA, and occasionally of exogenous origin. To address these challenges, we present exceRpt, the exRNA-processing toolkit of the NIH Extracellular RNA Communication Consortium (ERCC). exceRpt is structured as a cascade of filters and quantifications prioritized based on one's confidence in a given set of annotated RNAs. It generates quality control reports and abundance estimates for RNA biotypes. It is also capable of characterizing mappings to exogenous genomes, which, in turn, can be used to generate phylogenetic trees. exceRpt has been used to uniformly process all ∼3,500 exRNA-seq datasets in the public exRNA Atlas and is available from genboree.org and github.gersteinlab.org/exceRpt.

exRNA Atlas Analysis Reveals Distinct Extracellular RNA Cargo Types and Their Carriers Present across Human Biofluids.

To develop a map of cell-cell communication mediated by extracellular RNA (exRNA), the NIH Extracellular RNA Communication Consortium created the exRNA Atlas resource (https://exrna-atlas.org). The Atlas version 4P1 hosts 5,309 exRNA-seq and exRNA qPCR profiles from 19 studies and a suite of analysis and visualization tools. To analyze variation between profiles, we apply computational deconvolution. The analysis leads to a model with six exRNA cargo types (CT1, CT2, CT3A, CT3B, CT3C, CT4), each detectable in multiple biofluids (serum, plasma, CSF, saliva, urine). Five of the cargo types associate with known vesicular and non-vesicular (lipoprotein and ribonucleoprotein) exRNA carriers. To validate utility of this model, we re-analyze an exercise response study by deconvolution to identify physiologically relevant response pathways that were not detected previously. To enable wide application of this model, as part of the exRNA Atlas resource, we provide tools for deconvolution and analysis of user-provided case-control studies.

Integrative functional genomic analysis of human brain development and neuropsychiatric risks.

To broaden our understanding of human neurodevelopment, we profiled transcriptomic and epigenomic landscapes across brain regions and/or cell types for the entire span of prenatal and postnatal development. Integrative analysis revealed temporal, regional, sex, and cell type-specific dynamics. We observed a global transcriptomic cup-shaped pattern, characterized by a late fetal transition associated with sharply decreased regional differences and changes in cellular composition and maturation, followed by a reversal in childhood-adolescence, and accompanied by epigenomic reorganizations. Analysis of gene coexpression modules revealed relationships with epigenomic regulation and neurodevelopmental processes. Genes with genetic associations to brain-based traits and neuropsychiatric disorders (including MEF2C, SATB2, SOX5, TCF4, and TSHZ3) converged in a small number of modules and distinct cell types, revealing insights into neurodevelopment and the genomic basis of neuropsychiatric risks.